Bounds on gravitational wave backgrounds from large distance clock comparisons 
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Our spacetime is filled with gravitational wave backgrounds that constitute a fluctuating environ- 
ment created by astrophysical and cosmological sources. Bounds on these backgrounds are obtained 
from cosmological and astrophysical data but also by analysis of ranging and Doppler signals from 
distant spacecraft. We propose here a new way to set bounds on those backgrounds by performing 
clock comparisons between a ground clock and a remote spacecraft equipped with an ultra-stable 
clock, rather than only ranging to an onboard transponder. This technique can then be optimized 
as a function of the signal to be measured and the dominant noise sources, leading to significant 
improvements on present bounds in a promising frequency range where difi'erent theoretical models 
are competing. We illustrate our approach using the SAGAS project which aims to fiy an ultra 
stable optical clock in the outer solar system. 



I. INTRODUCTION 

The basic observables used for synchronizing remote 
clocks or ranging to distant events are built up on elec- 
tromagnetic signals exchanged between remote observers 
and compared to locally available atomic clocks |1| . This 
statement applies for example to planetary radar rang- 
ing ^ , lunar laser ranging , synchronizing orbiting 
clocks vifith Earth-bound standards Q or tracking and 
navigating probes in deep space @. 

Electromagnetic links feel the gravitational waves 
(GW) and this is currently the main route toward GW 
detection. It follows that GW affect ranging and Doppler 
tracking observables 0, [!]■ This effect has been thor- 
oughly studied in particular with the Pioneer and Cassini 
probes d, [l^, leading to constraints on the GW 
noise spectrum in some frequency range jT^ . These stud- 
ies constitute one of the windows on the physics of the 
stochastic GW backgrounds which permeate our spatio- 
temporal environment and have an astrophysical or cos- 
mological origin [l^ [l^, [H, [H, [13 ■ Their results have 
to be compared with bounds obtained through different 
observations [l^ (more discussions below). 

The aim of the present paper is to show that remote 
clock synchronization is also affected by GW and might 
be used to set new bounds. Timing is less sensitive than 
ranging at distances shorter than the GW wavelength, 
but this is no longer the case at large distances. Fur- 
thermore, the timing procedure can be arranged in order 
to get rid of uncertainties on the motion of the remote 
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clock, which might greatly improve the bounds on GW at 
low frequencies. The numbers will be discussed below by 
taking as an example the SAGAS project which aims at 
flying ultrastable optical atomic clocks in the outer solar 
system [l9j . These numbers heavily rely on the extremely 
good accuracy of modern atomic clocks [13, HU, ■ 

In the next section, we introduce and compare the ba- 
sic observables associated with ranging and timing. We 
then discuss their sensitivity to stochastic GW back- 
grounds as well as the noise sources involved in their 
measurement. We finally deduce the constraints on GW 
backgrounds which could be drawn from comparisons be- 
tween accurate clocks at large distances from each other 
in the solar system. 



II. RANGING AND TIMING OBSERVABLES 

We study the comparison between an atomic clock on 
board a probe and another one colocated with a station 
on Earth. The clocks are compared using up- as well as 
down-links. The uplink signal is emitted from ground at 
positions {ti, xi) in time and space and received in space 
at (i2,X2). The downlink signal is emitted from space 
at (t3,X3) and received on ground at (t4,X4) (see fig. [T]). 
The up-link is independent from the down-link i.e. ^3 — ^2 
can be chosen to take any value, positive or negative. 

The positions of the emission and reception events 
are connected by light cones, which may be calculated 
through a variety of theoretical methods (see for exam- 
ple psi. [2^. |2^] ) . These calculations give solutions which 
depend on the motions of Earth station and space probe, 
as well as on the gravitational field described by the met- 
ric. As the GW are weak modifications of the metric, we 
treat their effect as a perturbation of the solutions. 

These solutions may be written as relations between 
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FIG. 1: Principle of a general two-way link with t2 7^ ts 



clock indications corresponding to proper times elapsed 
on ground or space, say rf and r| for the uplink, and t| 
and t| for the downlink (see fig. [1]). We first define a 



ranging observable by 
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Td = tI 



(1) 



In the general case up- and downlinks are defined in- 
dependently from each other, but a specific configura- 
tion of interest is the situation where the links coincide 
at the space endpoint (t| = t| <^ ^2 = ^3)- Then 
represents a spatial distance corresponding to half the 
proper time elapsed on ground during the roundtrip of 
the signal to and from the probe. This special case is the 
only one realized in deep space probes so far, which were 
only equipped with a transponder that essentially only 
reflects the incoming signal. The range observable ([T|) is 
unaffected at lowest order by clock uncertainties in space 
(when using t2 = ^3) and can even be measured without 
a good clock on board. 

Another observable of great interest is the combination 
of Tu and Td defined with the opposite sign 



n 



(2) 



This timing observable ^ is unaffected at lowest order 
by uncertainties in the motion of the probe (when t2 = 
t^). For a probe equipped with a clock and a two-way 
system, one can choose to use either (or both) of the 
observables ^ and ^ with a free choice of the value of 
^3 — ^2 in order to optimize the measurement depending 
on the signal to be measured and the noise affecting the 
measurements. This is not the case for probes equipped 
only with a transponder, which are limited to the special 
case of ([T]) with t2 = H- 

We also introduce the time derivatives of ([T]) and ([2]) 



2/r 



yt = Tt 



-yd + Vu 



(3) 



The dot symbol represents here a derivation with respect 
to a commonly defined time t, chosen for any convenient 
argument. Note that is directly related to the Doppler 
tracking observable, which has been over the years the 
main source of information on the navigation of remote 
probes [1^. Meanwhile yt is directly related to the fre- 
quency comparison of distant clocks, the so-called syn- 
tonization observable. The variations of these quantities 
can be evaluated in the framework of a linearized approx- 
imation with a reasonably good approximation in the so- 
lar system. A more precise evaluation would be easy by 
using available methods [23j| and it would not change the 
qualitative discussions presented below. 

The variations in fSj are due to the effect of motion of 
the probe and ground station and perturbing effects like 
atmospheric delays, clock noise, etc... (see sect. 4) on 
one hand, and to the integrated effect of gravity along 
the propagation of the electromagnetic link on the other 
hand. The latter is given by an integral along the up- or 
downlink paths [u] and [d] 



STn,d = -TT 



u,d] 



dx' 



i.d ^^iiA 



da da 



da 



(4) 



hjj^ is the metric perturbation (/i^i, = g^j^^ — rj^i, with 
g^i, and rj^v the metric and Minkowski tensors) in the 
transverse traceless (TT) gauge; a is the affine parameter 
along the path measured as a length, and ^ the electro- 
magnetic wavevector reduced so that its time component 
is unity. 



III. SENSITIVITY TO STOCHASTIC GW 
BACKGROUNDS 

We now evaluate the effect of stochastic GW back- 
grounds as sources of noise on the electromagnetic links. 
For simplicity, we consider the background to be isotropic 
and unpolarized. We first recover well known results for 
the ranging case with t2 = is and then discuss the timing 
observable and the general case t2 7^ ts- 

To this aim, we introduce a plane wave decomposition 
of the GW background with wavevector k = fcn (fc = cj/c 
the modulus of k and n its direction) 



/i//(i,x) 



H- C.C. 



(5) 



We then write the Fourier component 5t[uj\ of the phase 
shifts appearing in eq.(|4]). At each frequency cj, this com- 
ponent is an integral over the direction n of GW plane 
waves with wavevectors k = {uj/c)r\ weighted by sensitiv- 
ity amplitudes. Considering as an example propagation 
along axis 1, we get 



5t[u] 



2nc^ 



d^nhjjjkn] 
An 1 - 



P[kn] 



(6) 
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The sensitivity amplitudes f3[kn] depend on the fre- 
quency (jj and the parameter /x = rii , which is the compo- 
nent of n along the direction of propagation of the elec- 
tromagnetic signal (here the axis 1). For a signal emitted 
at (ti, xi) and received at (^21X2) (uplink on figure 1) the 
sensitivity amplitude is given by [23| 

—i 

We will consider for simplicity the case of a stationary, 
unpolarized and isotropic background. The background 
may thus be characterized by a spectral density S'gw [w] 
giving the strain noise at a space point x (28j 

{hj^{t,^)hj^{0,^)) = -j^ _^5Q^[,^]cos(L^t) (8) 

The fluctuations of r are finally characterized by a noise 
spectrum Sr[uj] such that [2^ 

{dT{t)6T{0)) ^ f —Sr[to]cos{iut) 
Jo 27r 

5 

St[^] = ^ b[uj]SGwM (9) 

The dimensionless function b is obtained by averaging 
|/?[fcn]p over the direction n of the GW wavevector 

6H^(|/3MP)„ = f'^mn]\' (10) 
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FIG. 2: b (red-dotted) and c (blue-solid) as functions of a; = 

luT/tt 



For obvious symmetry reasons, it has the same form for 
up- and downlinks. 

The variation of b versus the dimensionless parameter 
X = loT/tt is shown as the red-dotted curve on Figure 2. 
The blue-solid curve on Figure 2 represents the covari- 
ance function c which describes the correlation of the up- 
and downlinks (* denotes complex conjugation) 

cH EE i (/3„[A;n]/3^[fcn] + /3d[fcn]/3:[A;n])„ (14) 
The average over n can be evaluated as 



A. Up- and downlinks with t2 = 

We first discuss the special case t2 = ^3- The sen- 
sitivity amplitudes for the up- and down- links are ob- 
tained directly from ([7]) by fixing the origin of coordinates 
t2 = ^3 = X2 = X3 = yielding 

^ l_t (1 ^ e-^(i+'')-^) (11) 

We have introduced shorthand notations for the propa- 
gation time T and the cosine /i (defined for the uplink) 

T = t2-ti=U-h (12) 

Throughout the paper we assume that the relative mo- 
tion of the space probe and the Earth during signal prop- 
agation is neghgible (so that ^2 — = ^4 — ^3)- More 
general expressions for the situation where this is not 
the case can also be obtained, but are beyond the scope 
of this work. The up- and downlink expressions are ex- 
changed by taking opposite signs for the cosine ^, the 
propagation time T and the global expression. 
The resulting 6— function is already known p?} 

, , , ^( A sin (2a;T) - 2wT\ 



cM = -47[cj]cos(wT) (15) 
, , cos(tjr) cos(tjT) sin(wr) 

We conclude this section by discussing qualitatively 
the shapes of the two curves b and c. We first notice that 
(3-a and /3d tend to become identical at the limit of short 
distances or low frequencies (/3u — /3d — (l ~ M^) "^"^ 
for LoT <C 1), so that b and c show the same behaviour 
~ ^cj^T^. The noise spectrum ^ is thus reduced to 
the simple form (corresponding to 5t^ ~ 5t^ ~ —-^hu) 

5,M = ^^GwH , ^T«l (16) 

At the high-frequency or large distance limit luT 3> 1 in 
contrast, b goes to a constant so that 

^-M = A^GwH , wT»l (17) 

Meanwhile, the correlation c between up- and downlinks 
remains sensitive to the distance even at large distances. 
This is simply due to the fact that (/3u/32 -I- /3d/3*) /2 con- 
tains a part (/i^ — l) (1 -t- cos {2ojT)) which is not blurred 
by the integration over /j. We also note that c, which is 
positive at low frequencies {ujT <C 1), is negative at high 
frequencies {ojT > f ). 
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B. Ranging and timing with t2 = 

We repeat now the same discussion in terms of the 
ranging and timing observables. It is clear from ([T]) and 
(12) that one can write expressions similar to ^ for St^ 
and Srt with the following sensitivity amplitudes 

A.^Jil* , A = zAi±A ,18, 

The noise spectra have the same form as © with the 
sensitivity functions 

KP] ^ {mkn]f) ^ , 6tM = (l/?t[fcn]l')^ (19) 

which can be written in terms of the already discussed 
functions b and c 

, - , b\uj] + cluj] , . , bluj] — clu] , , 

fcrM-^^Y^ , hM = ^ \ ^ ^ (20) 

It has also to be stressed that the correlation between 
the ranging and timing variables vanishes, as can be 
shown through an explicit calculation. As a matter of 
fact, the sensitivity amplitudes can be written as 

/3r = (sin cjT — /i sin /iwT 

-i^ (coscjT - cos [iijjT)) exp'''"^'^ 
(3t = (sin fj,ujT — fi sin loT 

+i {cos ujT - cos fj.ujT))exp'''"^'^ (21) 

and it turns out that {Pi (3^ + (3tP*) /2 is odd in /i and 
vanishes after the angular integration. Alternatively the 
fact that the correlation between St^ and (5rt vanishes can 
be directly inferred from the already discussed property 
b-a = ^d, which was attributed to a symmetry between 
up- and downlinks. It means that St^ and Srt appear as 
intrinsic and independent stochastic fluctuations of the 
positions in space-time of the end-points. 

The explicit expressions of the functions 6r and 6t can 
be obtained through an explicit calculation or alterna- 
tively deduced from (|13p and 



3 - cos {2ujT) 3 + cos {2ujT) 2 sin (2wT) 



5 + cos {2ujT) 1 



{uoTf 
cos {2u)T) 



{coTf 



{luTY 



(22) 



The associated plots are shown for b^ (blue-solid) and b^ 
(green-dashed) on Figure [3l 

Since /3u and /3d tend to become identical at the limit 
of short distances or low frequencies, it follows that the 
GW affect essentially the ranging observable. We effec- 
tively obtain in the limit wT 1 a much larger value 
for b, ~ 8(wr)Vl5 than for bt ~ 2(a;T)'^/15. This is 
obviously the reason why this case has been much more 
studied than the timing case. But is also clear that this 
is no longer the case for arbitrary frequencies. In partic- 
ular we know that c is negative at frequencies loT > -I so 
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FIG. 3: bi (blue-solid), 6t (green-dashed) and 6opt (red-dotted, 
see section 3.3) as functions of a; = loT/tt 



that bt is larger than b^- in this frequency range. We also 
notice that the oscillations in for and bt persist at large 
frequencies or large distances with simple behaviours 
b, ~ 1 - cos(2wr) /3 and h ~ 5/3 + cos (2a;T) /3 at 
ujT ^ 1. While this merely reflects the discussion al- 
ready devoted to c it is worth noting that the oscillations 
tend to disappear in the sum of b^ and bt that is also the 
single-link expression b. 



C. General case t2 ^ 

For a spacecraft equipped only with a transponder the 
only observable that can be obtained is the ranging 
(or its derivative j/r) for the special case t2 = (cf Fig. 
1). An onboard clock enables one to use also the timing 
observable rt, and additionally provides the possibility 
to freely choose the ground and onboard measurements, 
which are combined to form or rt, ie. to freely choose 
the value of t^ — t2- We will show below that this choice 
can be used to optimize the sensitivity of the observables 
to a particular signal (in the present case the GW back- 
ground) for a given measurement noise spectrum. 

The general expression ([7]) for the one-way sensitivity 
amplitude can be used directly in pS]) and p9|) , to obtain 
the sensitivity functions 6r,t[t^] for an up- and downlink 
separated by = t^ — t2 



b,[uj] = -- 27 cos (cjTsa + ujT) 

bt + 27 cos {ijT23 + LOT) (23) 



with b given in p3|) and 7 in p^ . The special case 
^2 = ^3 (equations (|22p) is recovered when setting = 
in eq. (|23p . These general sensitivity functions can now 
be used in the data analysis to choose the optimum value 
of T23 for each frequency w as a function of the link noise 
and signal travel time T. Figure 3 shows the optimized 
sensitivity function 6opt ^ calculated by maximising either 
of the equations One clearly obtains 6opt > b^.bt- 
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IV. MEASUREMENT NOISE 

In the following, we take the SAGAS (Search for 
Anomalous Gravity using Atomic Sensors) project 
as an example to illustrate the advantages and versatility 
provided by missions with an onboard clock and indepen- 
dent up- and downlinks. The SAGAS project proposes 
to fly a highly stable and accurate optical atomic clock 
and atomic accelerometer on an escape orbit in the so- 
lar system, up to a distance of 50 AU and beyond. It 
will use a continuous optical link for clock comparison, 
navigation, and data transfer, together with an X-band 
radio link as a backup. Science objectives are centered 
on tests of fundamental physics, in particular gravity on 
solar system scales and the exploration of the outer solar 
system, in particular the Kuiper belt. 

The optical link uses continuous transmission of a laser 
at /o ~ 444 THz in both directions (up and down) with 
1 W at emission, a 40 cm telescope on board the satellite 
and 1.5 m telescopes on the ground. Numerous pertur- 
bations on the link (atmospheric and instrumental losses, 
received photon flux and shot noise, stray light, etc..) are 
discussed in [l^, section 3.3.4. The fundamental science 
measurements of SAGAS are the frequency difference be- 
tween a local laser (optical clock) and an incoming laser 
beam at the same nominal frequency /o, both on board 
(up-link) and on the ground (down-link), sampled at 0.01 
Hz (see fig. [IJ. The measurements thus correspond to 
the observables tJu and j/d defined in Q. Including only 
terms whose noise contribution plays a significant role 
they can be expressed as 

_ fs{t2) - fsih) , Vs(^2)-Vg(tl) 

Jo c 
-I-Aytropo(ii) 

kiU) ~ Uih) , Vs(t3)-Vg(^4) 

yd = J 1- J-^d • 

Jo C 
-hAytropo(i4) (24) 

where /s,g are the frequencies of the space/ground 
laser (optical clock), Vg^g the associated velocity vec- 
tors, Nu,d the direction vectors of up- and down- 
links f N„ = , Nd = ;^°fe)-x,(t4) \ ^ 

V l|xB(t2)-Xg(ii)|| ' a ||xs(i3)-Xg(t4)|| y 

Aj/tropo the frequency change of the signal due to it cross- 
ing the Earth's troposphere. 

The noise coming from the different terms in ([24]) can 
be described equivalently by a power spectral density 
(PSD) Sy{f) or an Allan variance <Jy{T). In the following 
we consider the simple cases of a white frequency noise 
(terms proportional to /iq) and of a flicker frequency noise 
(terms proportional to /i-i) with the translation rule ^29l] 

crl{T) = ^ + 2\n2h^^ (25) 

Here the frequency / is in Hz and the integration time 
T in s (we have kept the notations and used in 



time and frequency metrology and which should not be 
confused with the metric perturbations hj^). We first 
discuss the different noise contributions on a single link 
and then evaluate the noise on the combined observables 
?/r and yt defined in ^ for the general case t2 ^3. 

A. One-way link measurement noise 

The space clock fractional frequency stability, as spec- 
ified in [l9|, corresponds to a white frequency noise 
(Tj,(r) = 1 X 10^ ^"'t^^/^ and an accuracy of iJy{T) = 
1x10^^^, ie. a flicker frequency noise at or below that 
level. The fractional frequency power spectral density 
(PSD) is then 

SyXI) - (2 X 10^28 + 7.2 X 10"3^/-i)/Hz (26) 

Though the flicker noise is likely to be lower than the 
projected accuracy, we use this conservative estimate. 

The ground clock stability is likely to be significantly 
better than the space clock stability by the time the 
mission is launched. Best present stabilities of opti- 
cal frequency standards are already below <Jy{T) = 3 x 
lQ-i5^-i/2 p2], with accuracies at 3 x lO^^"^ ;21]. Fur- 
ther rapid improvement of these numbers is expected. 
We therefore estimate the ground clock noise by the time 
of mission operation at ^^(r) = 5 x 10^ ^^r^^/'^ with a 
flicker component at 3 x 10"^®, so that the PSD is read 

SyM) = (5 X 10"'' + 6.5 X 10-36/-i)/Hz (27) 

The noise on the spacecraft velocity of SAGAS is de- 
termined by the integrated noise of the on-board ac- 
celerometer. Although orbit modeling is likely to im- 
prove on the raw accelerometer noise at low frequency, 
we use that as our conservative estimate for the purpose 
of this work. The accelerometer noise specified in [l^ 
is o-Q(r) 9 X 10~^°T^^/^ m/s^ per axis for 3D mea- 
surements and \/3 less when measuring only along the 
direction of signal transmission, of interest here. An ab- 
solute accelerometer based on cold atom technology is 
used to avoid long term drifts and biases (see [l^, sec- 
tion 3.1 for details). The expected absolute accuracy is 
5 X 10~^^ m/s^ taken again as the upper limit of the 
flicker acceleration noise. This translates into a velocity 
PSD of 

S^^/c = (1.5 X IQ-^'^f-'^ + 5.1 X lQ-^'^r^)/nz (28) 

For radio-frequency Doppler ranging, one of the domi- 
nant noise sources at low frequency is the uncorrected 
motion of the 34 m DSN antenna and of the station 
location itself [33|. For the optical link ground tele- 
scopes the motion of the mirror is likely to cause less 
of a problem, however the site movement plays a similar 
role as in the radio-frequency case. At high frequencies 
the motion of the ground station can be corrected us- 
ing gravity measurements, with best presently achieved 
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measurement noise levels of about 4 x 10~^^ m^/s^/Hz 
when using superconducting gravimeters [Slj. Alter- 
natively, positioning using global navigation systems 
(GNSS) and/or Satellite Laser Ranging (SLR) and/or 
Very Long Baseline Interfcrometry (VLBI) achieve sub- 
cm uncertainties [13] . Typically, GNSS positioning shows 
flicker noise in position over a wide range of frequen- 
cies [H, 113, HI], with best results at present at about 
Sxif) = 1 X 10~^/~^ m'^/Hz. The noise level shows some 
dependence on the number of visible satellites, and there- 
fore further improvement is expected with upcoming ad- 
ditional GNSS systems ^ . SLR and VLBI show white 
positioning noise 3j| , but at higher levels than GNSS at 
the frequencies of interest here (10^^ — 10^^ Hz). Com- 
bining local gravity measurements with GNSS position- 
ing we obtain as the minimum noise on ground velocity 




Frequency (Hz) 

FIG. 4: Power spectral densities of the dominant noise sources 
on a one-way link: Sy^ (magenta-crosses), Sy^ (light blue- 
dashed), Sy^/c (dark blue-solid), S^^/c (red-dotted), Sytropo 
(green-dash-dotted) . 



-Llx 10-36£^, />1.4x10-5Hz 
^ Hz 

S.,j, = 4.4 X 10-^^-^ , / < 1.4 X lO-^Hz (29) 

In (|29p we have assumed that superconducting gravime- 
ters display white noise down to 10~^ Hz, whereas the 
spectra shown in [3l| only show white noise down to 
10^'^ Hz, with the measurements being dominated by 
natural gravity fluctuations (the signal to be measured) 
at lower frequencies. This has to be considered as a pre- 
liminary estimation of low frequency noise unknown at 
present. 

Typically tropospheric delay models at optical frequen- 
cies have millimetric accuracy. Furthermore mapping 
functions that determine the variation of the delay with 
elevation (of interest here, the constant part of the de- 
lay playing no role on the frequency measurement) have 
been developed to sub-millimeter accuracy and success- 
fully tested on SLR data [s^j ■ Assuming that the residu- 
als from such models show white phase noise at the 1 mm 
level at 10 Hz sampling (typical pulse rate of SLR and 
LLR stations) we obtain a frequency PSD of 

Vopo(/) = 8.7x10-23/VHz (30) 

Note, that we pessimistically ignore correlations at high 
frequencies, which, given the slow motion of the satellite 
in the sky, should lead to decreased high frequency noise. 
Figure S] summarizes the noise sources on a one way link 
discussed above. 



B. Measurement noise on ranging and timing 

The noise on the ranging and timing observables de- 
fined in ([3]) is a combination of the noise affecting the in- 
dividual links, which are clearly correlated when forming 
the ranging or timing observables. Taking into account 



those correlations for the general case t2 7^ ^3 we obtain 

SyM) {(1 - cos(2^/Ti4)) Sy^if) (31) 

+ {l-cosi27r fT23))SyJf) 

+ (l+COs(2^/T23))5,,/c(/) 

+ (l + cos(2^/Ti4))5„jc(/) 
+ (l + cos(2^/Ti4))5,,,„,„(/)} 

SyM) ^\ {(1 + cos(2^/Ti4)) Sy,Sf) (32) 

+ (l+COs(2^/T23))5j,,(/) 
+ (1 - C0S(2^/T23)) S,^,,{f) 
+ (1 - cos(2^/Ti4)) S,j,{f) 
+ (l-cos(27r/Ti4))5,,„^„(/)} 

where we have defined Tij = tj — ti (cf. fig. [T]) . 

Figures 5 show the resulting noise PSD of the timing 
and ranging observables for the special cases T23 = 
(coincidence of up and down signals at the satellite) and 
2^14 = (coincidence of up and down signals on the 
ground). More generally, (|3ip and p2|) can be used to 
obtain the noise spectra of y-^ and yt for arbitrary values 
of T23 and ri4, with the constraint T14 — 123 = 2T. 

We note that the noise levels are significantly different, 
in particular at low frequency, illustrating the potential 
gain one can expect from using the optimal observable. 
In particular the observable y^ with T23 = used in " clas- 
sical" Doppler ranging shows several orders of magnitude 
larger noise at low frequency than the timing observable 
yt with T23 — 0. This is due to the cancelation of the 
onboard acceleromcter noise in (j32p leaving only the on- 
board clock as the dominant noise contribution at low 
frequency. This advantage has to be weighed against the 
different sensitivity functions &r and 6t as shown in fig. [3] 
in order to determine the optimal observable as a func- 
tion of Fourier frequency, which will be the subject of the 
following section. 
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FIG. 5: PSD of the noise affecting the two observables yt 
(upper graph) and j/r (lower graph) for the special cases ~ 
(blue-solid) and T14, = (red-dotted), and for a spacecraft 
to ground distance of 10 AU. 



V. CONSTRAINTS ON GW BACKGROUNDS 

The case of ranging or Doppler tracking (j/i ) has been 
discussed in numerous papers with the best bounds given 
by Doppler tracking of the Cassini probe [H, [13] ■ One- 
way hnking has also been studied since it is involved in 
the extremely impressive bounds derived from pulsar tim- 
ing [38, 39]. Here we will focus our attention on the case 
of synchronization between remote clocks. As made clear 
by the discussion of the preceding section, this points to 
experiments with excellent clocks at large distances, and 
we will take the SAGAS project as an example. We be- 
gin by discussing a somewhat simplified case illustrating 
the advantages of the different observables and combina- 
tions, and then go on to derive limits using the complete 
SAGAS noise sources as discussed in section 4. 

In order to discuss the attainable performances, we 
use the spectra associated with time derivatives of the 
phaseshifts induced by GW backgrounds (eq. [1]). These 
have to be compared to the phase variation induced by 
the noise sources discussed in section 4. The equality of 
the two provides the obtainable upper limit of S'gw 

5 

5"^ [uj] = uj'^Sr [uj] = -b [uj] Sow [1^] (33) 
o 

with Sy and b having different forms for the different 
observables (cf. equations (|23p . ((3T|) . ([32l)). 

For comparison with known bounds [l8|, we describe 
the gravitational noise in terms of the reduced grav- 
itational energy density f^cw commonly used to dis- 
cuss the cosmic backgrounds {Hq ~ 71 kms~^Mpc~^ ~ 
2.3 X 10^1** Hz is the Hubble constant) 



"GW — o u2 '-'GW 

6Hq 



(34) 



GW 



16^^ 
b 



(35) 



Collecting ([55)1 and , we deduce the expression of the 
bound obtainable on Qgw from that of the sensitivity 



A. Illustration using a simplified case 

To illustrate how the obtainable limits can be opti- 
mised using the available measurements and resulting ob- 
servables, we first consider an idealized case where only 
three noise sources play a significant role, the onboard 
clock and accelerometer, and the ground clock. Further- 
more, we will assume that all three noises consist of only 
white noise, at the levels indicated by equations (PSI [25] 
and[77|). ie. we will only consider the first terms of those 
equations in the expressions pTl [5^ for Sy , with all other 
terms set to zero. 

Figure 6 shows the resulting limits on fJcw as a func- 
tion of frequency for two satellite to ground distances: 
6 AU, the distance of Cassini when the GW experiment 
was carried out [30| . and 53 AU, the maximum distance 
envisaged for the SAGAS mission [l^ . Limits are shown 
for three observables: fir is obtained using the "classi- 
cal" Doppler ranging observable as defined in ^ with 
t2 = ^3- We recall that this is the only observable avail- 
able on space-probes equipped only with a transponder 
(the case of all probes flown so far) . ilt is obtained using 
the timing observable defined in ^ again with the con- 
dition t2 = ^3. r^opt is calculated by adjusting T23 in (|23[ 
1311 and I32p for each frequency in order to minimize the 
obtained limit on S^gw- 

In doing so, one can use either the ranging or timing 
combination, the obtained optimal limits being identi- 
cal (albeit for different values of T23). That property is 
the result of the periodic dependence of equations ([23l 
[3II and [32]) on T23 , which means that at any given fre- 
quency one can find two values of T23 for which the rang- 
ing and timing combinations yield the same limit on f^cw 
in equation ([35]) . However, the assumption that the up 
and down travel times are similar (see section 3.1) limits 
the allowed range of r23. Taking into account the maxi- 
mum relative probe-earth velocity (« 50 km/s) we limit 
r23 in the calculation of Oopt so that the up and down 
travel times do not differ by more than 1%. We then 
choose as fiopt the lower of the two limits obtained from 
j/r and yt with a free choice of T23 within the 1 % limit 
mentioned above. 

The graphs on figure[6]can be understood qualitatively 
and quantitatively when considering the expressions for 
the sensitivity functions ([23]) and the overall noise of the 
observables (PT|) and (|32p that enter into the calculation 
of r^GW in (|35| . We first discuss fir and fit in the case 
723 — 0, and then come to fiopt- 

At the low frequency hmit {coT « 1), 6r oc (luT)^ 
and bt oc (ujT)'^ (see section 3.2). The noise at low fre- 
quency is dominated by the space probe motion for Sy^ 
and is thus proportional to but that contribution is 
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FIG. 6: flaw versus /, for a probe to Earth distance of 6 AU 
(left) and 53 AU (right). The three curves show the hmits on 
Qr (blue-sohd), fit (green-dashed) and Qopt (red-dotted). 

entirely canceled in Sy^ because of the condition — 
in ([32]), so Sy^ oc cj^ (clock noise only). This leaves an 
overall lu~^ dependence for both fit and fl^, clearly dis- 
played in both plots of figure[n]at low frequency. Also, fly 
is significantly lower than fit at 6 AU, while the inverse 
is true at 53 AU (at low frequency). This is caused by 
a tradeoff between the difference in sensitivity functions 
and the involved noise sources. We find more specifically 
and to leading order 

lim ^ = 1 {u;Tf ^ (36) 

a;T<l fit A Sys 

where Sy^ and Sy^ are the first terms of ([^5]) and ([^5]) 
respectively, ie. Sy^/Syg ~ 3 x 10~^uj^^. Thus for large 
distances (greater than ~ 20 AU with SAGAS figures) 
the low frequency asymptote is lower for timing than for 
ranging, leading to the observed inversion of fl^ and fit 
when passing from 6 AU to 53 AU. 

At the other end of the spectrum {ujT >> 1), rang- 
ing outperforms timing (ie. fl^ < fit) at both distances. 
This can be easily understood when considering only the 
involved noise sources, as the sensitivity functions show 
oscillatory behavior and differ at most by a factor 3 (see 
section 3.2). At high frequency (between 10""* Hz and 
10~^ Hz), Sy^ is dominated by the space clock, but that 
contribution is entirely canceled in Sy^ because of the 



condition T23 = in (|3ip . leaving only a combination 
of space probe motion and ground clock noise. As at 
high frequency the space clock noise is significantly higher 
than that of the space probe motion or the ground clock 
(see figure \^ this leads to the observed advantage of fl^ 
over fit- Note also, the difference in slope between inl- 
and fit at high frequency, particularly visible at 53 AU, 
which can be easily understood from the slopes of the 
spectra of the different noise contributions (Sy^ cx J'^, 

Svs/c / ^). 

The lowest limits are obtained in the intermediate re- 
gion (luT ~ 1) with very different results for the two 
distances. As expected, even when using only the "clas- 
sical" ranging observable (f2r on figure [6]), limits improve 
with distance by about the ratio of distances (about an 
order of magnitude in the present case). However, it is 
clearly seen that almost another order of magnitude can 
be gained when taking advantage of the timing observ- 
able {fit in figure [5]), only available when using space 
probes equipped with an onboard clock and a two-way 
electromagnetic link. 

In that case, one can not only choose between y^- and 
yt, but also adjust the value of T23 in order to optimize 
the measurement for any given frequency. The result of 
such an optimization, Slopt, is shown in figure [6] As ex- 
pected it is below fl^ and fit at all frequencies and for 
both distances. Although the overall improvement is not 
spectacular, one obtains the "best of both worlds", in 
particular at 53 AU where flopt follows fit at low fre- 
quency and fly at high frequency. We notice a slight im- 
provement on fit at low frequency and 53 AU, which can 
be understood by considering the series expansion of the 
sensitivity function bt in (j23p . Additionally to the term 
in (wT)^ present in the case T23 — (see section 3.2) one 
now obtains a term proportional to (ujT)'^ (10X23)'^ which 
can be significantly larger. However, when T23 ^ 0, low 
frequency noise from the space probe motion is added, 
the trade off between the two leading to the small im- 
provement of flopt over fit seen on figure [S] 

In conclusion, the simplified case used in this section 
illustrates the advantages of having an onboard clock and 
a two-way link, which allows one to "fine tune" the data 
analysis as a function of the expected signal and the noise 
sources affecting the raw measurements. In this exam- 
ple, the sensitivity to GW backgrounds at large distance 
(53 AU) is improved by about an order of magnitude 
over the " classical" case fir by choosing the optimal com- 
bination of the available measurements on ground and 
onboard the space probe. Similar (up to a factor 20) 
improvements are observed when taking into account all 
noise sources discussed in section 4. 



B. SAGAS limits on GW backgrounds 

We now repeat the calculations described in the pre- 
vious section, using the example of SAGAS including 
all noise sources described in section 4. For clarity, we 
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show only the resuhing optimal limits flopt for a range 
of distances (see figure [7]). The Big Bang Nucleosynthe- 
sis (BBN) bound, which corresponds to a flat floor with 
f^BBN ~ 1-5 X 10^^ has been drawn for comparison. 




1Q-'L : : ^ : ■ 

10 10" 10" 10 10 



Frequency (Hz) 

FIG. 7: f2opt versus /, for probe to Earth distances (top to 
bottom) of 20 AU (light blue), 30 AU (magenta), 53 AU 
(green), 100 AU (red) and 200 AU (dark blue). The hori- 
zontal black line indicates the BBN bound. 

As expected, the sensitivity to GW is improved by go- 
ing to large distances. The frequency at which the lowest 
limits are reached is typically in the region where fT 1 
ie. decreasing with increasing distance. For SAGAS (pro- 
jected distance ~ 53 AU) the lowest obtained limit is 
about flow < 7 X 10~^ around 2 x 10^^ Hz and es- 
sentially determined from the timing observable. It is 
more than 3 orders of magnitude below the best directly 
measured limits in the 10^^ to 10^'' Hz band, obtained 
from the Cassini probe [s^, and about a factor 2 below 
the BBN bound. At larger distances the timing measure- 
ment could even approach the pulsar bound around a few 
10~®, but at significantly hig her frequencies (the pulsar 
bounds are at a few nHz |38|). 



C. Discussion 

Existing bounds on stochastic GW backgrounds are 
spread over a huge frequency range from 10~^^ Hz to 
10'^ Hz, corresponding to cosmological bounds obtained 
from measurements of the 3 K microwave background 
(COBE) at the lower end, and modern ground based GW 
detectors (LIGO, VIRGO) at the upper end (see e.g. fig- 
ure 14 of [IBl)- This large frequency range is patchily 
covered, with COBE limits at 10" - lO^^^ Hz, pul- 
sar bounds around 10^^ — 10^^ Hz, spacecraft Doppler 
ranging covering three orders of magnitude (10~^ — 10^"^ 
Hz), and GW detectors setting limits around 10^ Hz. 
This is complemented by an indirect upper limit derived 
from models of BBN which corresponds to a flat floor of 
^GW < 1-5 X 10-5 at all frequencies > 10-^*' Hz. The 
latter is already outperformed by the pulsar limits at low 
frequency and is expected to be outperformed at high fre- 



quency by LIGO and VIRGO measurements in the near 
future. 

In this landscape, limits obtained from spacecraft 
tracking play an important role as they fill a large part 
of the gap in frequency between the pulsar limits and 
those obtained from ground based detectors. Unfor- 
tunately, the obtained bounds are presently limited to 
^^GW < 0.025 30]. Improvements in this band will be 
particularly useful, especially when they will approach 
or surpass the BBN bound, as would be the case with 
future missions like SAGAS. As shown above, we expect 
that such missions could provide limits on Oqw down to 
parts in 10^® for SAGAS and below for missions at even 
larger distances. 

More generally, it is important to obtain experimen- 
tal constraints on ficw sA, all frequencies, as many of 
the models that predict such GW backgrounds are fre- 
quency dependent (eg. cosmic strings models, pre Big 
Bang models,...) and only poorly constrained in the 
presence of frequency gaps. In that respect the future 
space interferometric GW detector (LISA) plays an im- 
portant role, as it should provide extremely low limits 
(down to r^GW ^ 10"^^) in the still largely unconstrained 
frequency range of 10""' Hz to 10"^ Hz. It should even be 
able to observe the astrophysical GW background from 
an ensemble of galactic binary stars, estimated to be too 
low for any other present or planned detector, but within 
the reach of LISA. 



VI. CONCLUSION 

Doppler ranging to distant space probes provides the 
presently most stringent upper bounds on GW between 
10~^ and 10"'^ Hz. Those bounds are obtained by "pas- 
sive" ranging, where the space probe only serves as a 
"reflector" of the signal emitted from the ground. We 
have shown that the sensitivity can be significantly im- 
proved when having a clock onboard, so that the up and 
down signals are independent (asynchronous link) and 
can be combined in an optimal manner adapted to the 
signal to be measured and the noise affecting the link. 
We have derived explicit expressions for the sensitivity 
of all possible link combinations to a GW background. 
Using the example of the SAGAS project, we have evalu- 
ated the sensitivity of such a mission to GW backgrounds 
for optimal signal combinations and as a function of dis- 
tance, with a potential improvement by over three orders 
of magnitude on best present limits. 

Let us notice the similarities between the calculations 
of the present paper and those previously devoted to the 
effect of stochastic GW backgrounds on inertial sensors 
built on atomic interferometry [28, 40]. The sensors of 
interest in the present paper are the atomic clocks the 
indications of which are compared through electromag- 
netic links. As these links cannot be protected against 
the action of GW backgrounds, there exists an ultimate 
noise in clock synchronisation due to the presence of this 
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universal fluctuating environment. It has been shown in 
the present paper that timing can be more sensitive to 
this environment than ranging, provided that extremely 
large distances are considered, as it is the case in the 
SAGAS project. 

In our estimations we have chosen a conservative ap- 
proach where the noise on the spacecraft motion is de- 
termined solely by the measurement noise of the onboard 
accelerometer. Previous deep space probes, in particular 
the Cassini mission jsoj , did not have an accelerometer on 
board, and all non-gravitational accelerations acting on 
the probe where determined by fitting acceleration mod- 
els to the ranging data. The PSD of the residuals is most 
likely dominated by ground station and antenna motion 
at low frequency, and in particular around the diurnal 
frequency and its harmonics (see figure 1 of ,30]). For 
a mission like SAGAS this suggests an analysis strategy 
based on the cancellation of the ground station motion 
rather than that of the space probe. In frequency re- 
gions where ground station noise is dominant one would 
use the timing observable ^ giving rise to Sy^ (/) of ([5^ . 
but with the condition Tn ~ (coincidence of up and 
down signals at the ground antenna). As can be easily 
seen from (|32p . this leads to cancellation of noise from 
the ground station motion and the troposphere, leaving 
space clock instabilities and space probe motion as the 
dominant noise sources (see figure To evaluate the 
limits obtained in this scenario requires a more detailed 
investigation of the space probe motion, the effect of fit- 
ting acceleration models, the improvements in the fits 



from in situ acceleration measurements, etc., which are 
beyond the scope of this paper. Nonetheless, this alter- 
native approach well illustrates the versatility of using 
an asynchronous link that allows choosing the optimal 
data combination strategy, even after launch, and as a 
function of the observed noise levels. 

Finally, we point out that this data combination strat- 
egy can be adapted and optimized for any signal that is 
to be measured. The GW backgrounds discussed in this 
paper give rise to sensitivity functions (|23p which enter 
the parameter flow in ([55]) together with the noise (PT|) . 
P2p . That parameter is then optimized over a broad 
frequency range by varying T23 or T14. A similar proce- 
dure can be used for other science objectives (e.g. test 
of the gravitational time delay during occultation, mea- 
surements of planetary gravity, trajectory determination 
during fly by, etc..) by deriving the appropriate sen- 
sitivity functions and calculating the parameters to be 
optimized. It is likely to allow significant improvements 
in those measurements as well. 
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